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AN OPTICAL ENCRYPTION AND DECRYPTION METHOD AND SYSTEM 
FIELD OF THE INVENTION 

5 The invention relates to securing of information utilising imaging technologies and 
more specifically to phase and amplitude encryption and decryption of images. 

BACKGROUND OF THE INVENTION 

10 It is well known to form an image by phase contrast imaging methods in which phase 
modulation of light is converted into intensity modulation. As opposed to intensity 
modulation, phase modulation does not involve loss of energy. 

In published patent application no. WO 96/34307, which is hereby incorporated by 
15 reference, a phase contrast imaging method is disclosed for calculating phasor values 
of a phase mask for synthesising a desired intensity pattern. 

In "Securing Information with Optical Technologies", Bahram Javidi, Physics Today 
Vol. 50 f No. 3, March 1997, pp 27-32, a method and system for optically securing 

20 information is proposed. An image is encrypted in a 4f-lens optical configuration (i.e. 
comprising two Fourier transforming lenses) by Fourier transforming the image and an 
input phase mask with a first lens of the 4Mens system. A Fourier plane phase mask 
having phasor values e ib,a * P) is positioned in the Fourier plane of the first lens and an 
encrypted image is formed with a second lens that Fourier transforms images in the 

25 Fourier plane of the first lens. The encrypted image is decrypted in a similar 4f-lens 
configuration in which another Fourier plane phase mask (the key) having phasor 
values e" ^b,a • p, is positioned in the Fourier plane of the first lens. 

It is a disadvantage of the known cryptographic method that encrypting an image both 
30 in the object plane and the Fourier plane leads to generation of speckle patterns in the 
decrypted image thereby corrupting data having been encrypted. 
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It is another disadvantage of the known cryptographic method that extremely accurate 
three-dimensional positioning of the phase mask in the Fourier plane is required for 
successful encryption and decryption. 

5 It is yet another disadvantage of the known cryptographic method that both amplitude 
and phase have to be recorded in the encrypted mask. 

SUMMARY OF THE INVENTION 

10 It is an object of the present invention to provide an apparatus of the above kind 
which apparatus is robust, compact, and simple to design and relatively cheap to 
manufacture. 

It is another object of the present invention to provide an improved method and 
1 5 apparatus without the above-mentioned disadvantages. 

It is still another object of the present invention to provide an improved method and 
apparatus for optically securing data utilising phase contrast imaging. 

20 According to the invention a method is provided of decryption of an encrypted image 
having a non-encrypted image intensity pattern l(x\y'). 

The image is encrypted into a mask having a plurality of mask resolution elements 
(x m ,y m ) by encoding the image into the mask with an encoded phase value <t>(x m ,y m ) 

25 and an encoded amplitude value a(x m ,y m ), and by encrypting the mask by addition of 
an encrypted phase value (Mx m ,yJ to the encoded phase values <fr(x m( y m ) and by 
multiplication of an encrypted amplitude value a c (x m ,yj with the encoded amplitude 
value a(x m/ yj. Thus, each mask resolution element (x m ,y m ) modulates the phase and 
the amplitude of electromagnetic radiation incident upon it with the complex value 

30 a{x m ,yJa c ^x m ,y m )e^ m)+i ♦ c,,lm ^ , . 



The method comprises the steps of decrypting the mask by radiating electromagnetic 
radiation towards the mask, and inserting into the path of the electromagnetic 
radiation a complex spatial electromagnetic radiation modulator comprising modulator 
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resolution elements (x d ,y d ), each modulator resolution element (x d ,y d ) modulating the 
phase and the amplitude of electromagnetic radiation incident upon it with a 
predetermined complex value a d (x d ,Y d )e f ♦ d,, ' d • v^ ^ the decrypting phase value <Mx d ,y d ) and 
the decrypting amplitude value a d {x d ,y d ), respectively, of a modulator resolution 
5 element (x d ,y d ) being substantially equal to -<|> c (x m ,Y m ) and a c mt ix m ,y m ), respectively, of a 
corresponding mask resolution element (x m ,yj, and imaging the mask and the 
electromagnetic radiation modulator onto the image having the image intensity pattern 
l(x\y*). Throughout the present description the combination of the mask and the 
electromagnetic radiation modulator are denoted an encoder. 

10 

The complex spatial electromagnetic radiation modulator may be positioned anywhere 
in the path of the electromagnetic radiation. When the complex modulator is 
positioned adjacent to the mask, the encrypted mask is decrypted when the 
decrypting phase value <|> d (x d ,y d ) and the decrypting amplitude value a d {x d ,y d ), 

1 5 respectively, of a modulator resolution element (x d ,y d ) being substantially equal to - 
<t>c(Xm,Vm) and ac ^x^yj, respectively, of a corresponding adjacent and aligned mask 
resolution element (x m ,y m ). The complex modulator may be physically positioned 
remotely from the mask and imaged onto the position of the mask whereby the 
complex modulator is virtually positioned adjacent to the mask. The complex 

20 modulator may be positioned virtually or physically in another position than adjacent 
to the mask still performing the decrypting function when the amplitude and phase 
values of the complex modulator are transformed accordingly, e.g. by a Fresnel 
propagator. Thus, the mask and the decrypting complex modulator is aligned when 
corresponding resolution elements of the mask and the complex modulator are imaged 

25 onto the same desired resolution element of the image having the image intensity 
pattern l(x\ y'). 

It is an important advantage of the present invention that the encrypted image may 
selectively require recording in the mask of phase values or amplitude values or a 
30 combination thereof. 
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It is another advantage of the present invention that decryption is performed in a plane 
adjacent to the mask or an equivalent plane whereby generation of speckles in the 
decrypted image is suppressed. 

5 It is yet another advantage of the present invention that there is no requirement of 
positioning a radiation modulator in the Fourier plane whereby an accurate three- 
dimensional positioning requirement is avoided. 

It is still another advantage of the present invention that only a single key for 
10 encryption is required. 

It is preferred to encode the non-encrypted image into a phase mask and encrypt the 
phase mask by adding encrypting phase values to the encoded phase values. An 
encrypted phase mask is extremely difficult - if not impossible - to replicate by 
1 5 counterfeiters. Further, phase masks may be readily produced while masks requiring 
recording of both amplitude and phase are extremely complicated to produce, typically 
requiring production of two masks to be accurately superpositioned. The same applies 
to the decrypting complex modulator. 

20 The imaging may be performed with a common path interferometer, such as a phase 
contrast imaging system, a dark field imaging system, a field absorption imaging 
system, a point diffraction imaging system, a Smartt interferometer, a Schlieren 
interferometer, etc, and any combination hereof. 

25 According to a preferred embodiment of the invention a method is provided for 
reconstructing an intensity pattern l(x\ y*) , comprising 

Spatial modulation of electromagnetic radiation with the encoder for modulation of the 
phase of the incident electromagnetic radiation by phasor values of individual 
30 resolution elements of the encoder, each phasor value e i * U v) being determined in such 
a way that 

1) the values of the Fourier transformed phasors attains predetermined values for 
predetermined spatial frequencies, and 
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2) 



the phasor value of a specific resolution element of the encoder corresponds to 
a distinct intensity level of the image of the resolution element in the intensity 



pattern, 



5 Spatial phase filtering of electromagnetic radiation with a spatial phase filter for phase 
shifting of a part of the electromagnetic radiation, and 

Generating the intensity pattern l(x',y') by interference in the image plane of an 
imaging system between the part of the electromagnetic radiation that has been phase 
10 shifted by the phase filter and the remaining part of the electromagnetic radiation. 

Although, the present method is related to encoding in two spatial dimensions (planar 
encoding), the principles of the method may be utilised for phase encoding in one to 
three spatial dimensions and/or in the temporal dimension. 

15 

The electromagnetic radiation may be of any frequency range of the electromagnetic 
spectrum, i.e. the gamma frequency range, the ultraviolet range, the visible range, the 
infrared range, the far infrared range, the X-ray range, the microwave range, the HF 
(high frequency) range, etc. The present method is also applicable to particle radiation, 
20 such as electron radiation, neutron radiation, etc. 

Preferably, the electromagnetic radiation is monochromatic or quasi-monochromatic so 
that the energy of the electromagnetic radiation is concentrated in a narrow frequency 
bandwidth. Since the intensity pattern is reconstructed by interference of two 

25 electromagnetic waves emitted from a common source of electromagnetic radiation, 
the phases of which have been changed differently, it is required that the frequency 
range of the emitted electromagnetic radiation is sufficiently narrow to ensure that the 
two waves of electromagnetic radiation are coherent so that their superposition 
generates the desired intensity pattern. If the frequency range is too broad, the two 

30 waves will be incoherent and the phase information will be lost as superposition of 
non-coherent waves results in a summation of the intensities of the two waves. It is 
required that the difference between individual delays of electromagnetic radiation to 
be superpositioned is less than the wavelength of the radiation. This is a relaxed 
requirement that allows the electromagnetic radiation to be relatively broad-banded. 
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For example in the visible range a Xe-lamp or a Hg-lamp can be used as a light source 
in a system according to the present invention with the advantage compared to a laser 
light source that speckle noise is reduced. The requirements of the spatial coherence 
of the electromagnetic radiation depend upon the space bandwidth product of the 
5 corresponding system and how close the required system performance is to the 
theoretically obtainable performance of the system. 

Preferably, the electromagnetic radiation is generated by a coherent source of 
electromagnetic radiation, such as a laser, a maser. a phaselocked laser diode array, 
10 etc. However a high pressure arc lamp, such as a Hg lamp, a Xe lamp, etc, may also 
be used and even an incandescent lamp may be used as a source of electromagnetic 
radiation in a low performance system. 

The encoder changes the phase of an electromagnetic wave incident upon it. The 
1 5 encoder may transmit or reflect the incident electromagnetic wave. The encoder is 
divided into a number of resolution elements, each of which modulates the incident 
electromagnetic wave by changing its phase by a specific predetermined value. The 
predetermined values are assigned to each resolution element in different ways 
depending upon the technology applied in the component. For example in spatial light 
20 modulators, each resolution element may be addressed either optically or electrically. 
The electrical addressing technique resembles the addressing technique of solid-state 
memories in that each resolution element can be addressed through electronic circuitry 
to receive a control signal corresponding to the phase change to be generated by the 
addressed resolution element. The optical addressing technique addresses each 
25 resolution element by pointing a light beam on it, the intensity of the light beam 
corresponding to the phase change to be generated by the resolution element 
illuminated by the light beam. 

Spatial phase modulation may be realised utilising fixed phase masks, devices 
30 comprising liquid crystals and being based on liquid crystal display technology, 
dynamic mirror devices, digital micro-mirror arrays, deformable mirror devices, 
membrane spatial light modulators, laser diode arrays (integrated light source and 
phase modulator), smart pixel arrays, etc. 



A spatial phase filter is typically a fixed phase mask, such as an optically flat glass 
plate coated with a dielectric layer at specific positions of the glass plate. However, 
the spatial phase modulators mentioned in the previous section may also be used for 
spatial phase filters. 

5 

The imaging system maps the phase modulating resolution elements of the encoder on 
the target surface of the reconstructed intensity pattern. It may comprise a 4f-lens 
configuration (two Fourier transforming lenses utilising transmission of light or one 
Fourier transforming lens utilising reflection of light) or a single imaging lens. However, 
10 any optical imaging system providing a filtering plane for the spatial phase filter may 
be applied in a phase contrast imaging system. 

In the method according to the present invention, the reconstructed intensity pattern 
is generated by superposition of two electromagnetic waves in the image plane of the 

1 5 imaging system. The encoder changes the phase values of an electromagnetic wave 
incident upon it and the imaging system directs the electromagnetic wave with 
changed phases reflected from or transmitted through the encoder towards the spatial 
phase filter. The phase filter phase shifts a part of the electromagnetic radiation and 
the imaging system is adapted to superimpose in the image plane the phase shifted 

20 part of the electromagnetic radiation with the part of the electromagnetic radiation 
that is not phase shifted by the spatial phase filter. 

According to a preferred embodiment of the invention, the encoder is positioned at the 

front focal plane of a lens while the spatial phase filter is positioned in the back focal 
25 plane of the lens, whereby a first electromagnetic field at the encoder is Fourier 

transformed by the lens into a second electromagnetic field at the phase filter. Thus. 

specific spatial frequencies of the first electromagnetic field will be transmitted 

through the spatial phase filter at specific positions of the phase filter. For instance. 

the energy of the electromagnetic radiation at zero frequency (DC) is transmitted 
30 through the phase filter at the intersecting point of the Fourier plane and the optical 

axis of the lens also denoted the zero-order diffraction region. 

It is presently preferred that the spatial phase filter is adapted to phase shift the DC- 
part of the electromagnetic radiation and to leave the remaining part of the 
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electromagnetic radiation unchanged or, alternatively, to leave the DC-part of the 
electromagnetic radiation unchanged and to phase shift the remaining part of the 
electromagnetic radiation. The last alternative is preferred when the energy level of 
the DC-part of the electromagnetic radiation is so high that the phase shifting part of 
5 the phase filter will be destroyed by it. For example in laser cutting, the DC level of 
the laser beam can be so high that a phase shifting dot positioned at the intersecting 
point of the DC part of the laser beam at the phase filter would evaporate. It is also 
possible to block the electromagnetic radiation (no transmittance) in the zero-order 
diffraction region, however, the DC energy of the radiation is then lost. 
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Below, an expression of the intensity of the reconstructed intensity pattern as a 
function of the phasor values e*'*-* 1 of the encoder, when the DC-part of the 
electromagnetic radiation is phase shifted, is deduced. 

15 Electromagnetic radiation incident on the encoder can be described by a function 
A(x,y), where A(x,y) is a complex number (amplitude and phase) of the incident field 
on the point (x,y) of the encoder. At the point (x.y), the encoder modulates the phase 
of the incident radiation with a value <t>(x,y) so that the field after reflection by or 
transmission through the encoder may be described by the function A(x,y) # e M "\ 

20 e™*-" being the phasor value of the point (x.y) of the encoder. As A(x,y) preferably is a 
constant value over the entire surface of the encoder, the term is left out of the 
following equations for simplicity. 

The expression of the electromagnetic radiation incident on the spatial phase filter 
25 may now be separated into an AC-term and a DC-term. If the DC-term of the field is 
denoted a , the AC-term of the field is given by the term - a . Since the spatial 

phase filter changes the phase of the DC-part of the electromagnetic radiation by 6, 
the intensity of the reconstructed intensity pattern at the image plane of the imaging 
system is given by: 

30 

l(x\y') = I a*"'" 1 + a(e ,8 -1)l 2 
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wherein (x\y*) is the co-ordinates of the image of the point (x,y) of the encoder 
formed by the imaging system in the image plane. 



It should be noted that the second term of the equation is a complex number that 
5 adds to the phasors e*'* ¥ ' of the encoder and may be interpreted as a contrast control 
parameter for the reconstructed intensity pattern l(x',y'). 

According to a preferred embodiment of the invention, the average value of the 
phasors is adjusted in order to control the range of intensity levels. 

10 

Instead of phase shifting the DC-part of the electromagnetic radiation, it is also 
possible to synthesise a prescribed intensity pattern by phase shifting other parts of 
the electromagnetic radiation by adapting the spatial phase filter to phase shift 
electromagnetic radiation incident upon one or more arbitrary regions of the phase 
1 5 filter and leaving the phase of the remaining part of the electromagnetic radiation 
unchanged and then superimposing the two parts of the electromagnetic radiation. 
The corresponding mathematics and the corresponding design procedures for the 
encoder and spatial phase filter will of course be more complicated than for the 
method described in the previous section. 

20 

A simple example of phase shifting a part of the electromagnetic radiation of a spatial 
frequency different from the zero frequency is provided by moving the DC-part of the 
electromagnetic radiation to another spatial frequency in the Fourier plane (identical to 
the plane of the spatial phase filter) utilising an optical component with an appropriate 
25 carrier frequency (i.e. a grating or a prism) or, preferably, encoding the function of a 
grating or a prism into the encoder, and adapting the spatial phase filter to change the 
phase of the electromagnetic radiation at this spatial frequency and to leave the phase 
of the remaining part of the electromagnetic radiation unchanged. 

30 According to another preferred embodiment of the invention, the encoder is not 
positioned in the back focal plane of the lens but in the Fresnel region of the lens 
instead. In this case, the electromagnetic field at the phase filter will be given by a 
Fresnel transformation of the electromagnetic field at the encoder. This further 
complicates the mathematics and the design procedures, for example the term in 
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equation (7) has to be substituted by the value of the Fresnel transformation at the 
point(s) of phase changes of the phase filter. However, the Fresnel transformation 
may be calculated from a Fourier transformation by multiplication of the phasor values 
of the encoder by a quadratic phase factor followed by a Fourier transformation. 

5 

It is an important aspect of the present invention that each intensity level of the 
reconstructed intensity pattern for each resolution element may be generated by at 
least two different phasor values of a resolution element of the encoder. 

10 For example, when the spatial phase filter phase shifts the DC-part of the 

electromagnetic radiation, it will be shown later that, advantageously, the average a 
of the phasors of the resolution elements of the encoder should be equal to Vz and the 
value of the phase shift 9 should be equal to n In this case, the intensity of the 
reconstructed image pattern at the image (x\y') of the resolution element (x,y) will be 

1 5 given by: 

l(x', y') = 2(1 - cos 4>{x\y , H 

It is seen that complex conjugate phasors (values of $ of opposite sign) result in 
20 identical intensity levels l(x\y'). It can be shown that for any value of the modulus of 
the average of the phasors |a| two phasors exist that will generate identical intensity 
levels of the reconstructed intensity pattern. 

Further, if the spatial phase filter phase shifts parts of the electromagnetic radiation 
25 different from the DC-part, the phasor value that generates a specific intensity level 
will depend on the position of the resolution element in question, i.e. the phasor value 
and the position of the resolution element with that phasor value together define the 
intensity level at the image of the resolution element in the reconstructed intensity 
pattern. Still, it is true that for each resolution element of the encoder, each intensity 
30 level of the reconstructed intensity pattern may be represented by one of two 
different phasors of complementary phase values. 
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This freedom of being able to select, for each intensity level to be generated and for 
each resolution element of the encoder, one of two phasors is used to control the 
phase of the Fourier transform of the phasors at specific spatial frequencies by 
selection of phasors with appropriate phase values to ensure two intervals of biunique 
5 functional dependence between phasor values and corresponding intensity values. 

This freedom of choice of phasors may be utilised to select phasors of neighbouring 
resolution elements of the encoder with a maximum difference between them, thereby 
generating an electromagnetic radiation emitted from the encoder with a maximum 
1 0 content of high spatial frequencies which will generate a good separation of the DC 
part of the electromagnetic radiation from its AC part. However, any other strategy of 
selecting between two possible phasor values of each resolution element may be 
chosen to generate a desired spatial frequency content of the electromagnetic 
radiation. 



15 



20 



Preferably, the phase of the Fourier transform of the phasors at specific spatial 
frequencies is adjusted in order to control whether the relation between each phasor 
and the corresponding intensity level is a monotonic increasing or a monotonic 
decreasing function. 

Below, a set of different methods are described that are provided according to the 
present invention for adjustment of the modulus of the Fourier transform of the 
phasors at specific spatial frequencies to attain a prescribed value. If convenient, the 
methods may be combined. 
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According to one of the methods, the individual phasors of the resolution elements of 
the encoder are adjusted by a constant value until the desired value of the modulus of 
the Fourier transform of the phasors at specific spatial frequencies is attained while 
maintaining prescribed relative intensity levels between intensities of resolution 
30 elements of the intensity pattern, i.e. iteratively. 

According to another method, the individual phasors of the resolution elements of the 
encoder are adjusted utilising histogram techniques known from image processing. A 
histogram is a bar chart showing the number of resolution elements of the 
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reconstructed intensity pattern with a specific intensity value as a function of the 
intensity value. Any histogram technique, such as histogram equalisation, adapting the 
histogram to a predetermined distribution, etc., may be used iteratively until the 
modulus of the Fourier transform of the phasors at specific spatial frequencies attain 
5 the prescribed value. 

According to yet another method, the phasor pattern of the encoder is spatially scaled 
in order to adjust the modulus of the Fourier transform of the phasors at specific 
spatial frequencies. 
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According to still another method, the modulus of the Fourier transform of the phasors 
at specific spatial frequencies is adjusted utilising half tone coding techniques, such as 
raster techniques, area ratio modulation, spot diameter modulation, etc. 

15 It is seen from the description above that the intensity levels may differ from one 
reconstructed intensity pattern to the next as a consequence of the adjustments of 
the modulus of the Fourier transform of the phasors at specific spatial frequencies. 
Thus, it is preferred to control the power of the radiation source in dependence of the 
intensity range of the intensity pattern so that a sequence of different intensity 

20 patterns shows uniform intensity levels. 

According to a preferred embodiment of the invention, the shape of the phase filter is 
adapted to match the spatial frequency content of the phasors of the encoder, e.g. to 
optimise the desired separation of the part of the electromagnetic radiation to be 
25 phase filtered from the remaining part of the electromagnetic radiation. 

It is within the scope of the present invention that the imaging system further 
comprises zooming means for variable scaling of the reconstructed intensity pattern. 
The zooming of the imaging system may be dynamically controllable, e.g. in response 
30 to the scaling of the pattern of phasor values of the encoder. 

According to the present invention, the power of the radiation source may be 
controllable in response to the spatial scaling of the pattern in the encoder and/or the 
zooming of the focusing system. 
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In order to provide a compact and integrated system according to the present 
invention, the optical function of a Fourier-transforming lens is encoded into the 
phasors of the encoder. The Fourier transforming lens may be refractively or 
diffractively encoded. 

5 

Similarly, the optical function of an output lens may be encoded into the phase f.lter 
either refractively or diffractively. 

Further, compensation may be encoded into the phasor values of the encoder so that 
10 part of the electromagnetic radiation modulated by the encoder has a substantially flat 
intensity profile in the image plane. Without this compensation, part of the 
electromagnetic radiation modulated by the encoder will have a flat profile with 
perturbations resulting from the phase filtering superpositioned upon it. This may 
cause "ringings" (oscillations) at the edges of the reconstructed intensity pattern. 

15 

According to another preferred embodiment of the invention, the source of 
electromagnetic radiation comprises one or more light sources of different 
wavelengths corresponding to three different colours, such as red, green and blue, for 
generation of intensity patterns of arbitrary colours. Further, several independent 
20 systems each one illuminated by its own wavelength can be combined into a single 
multi-wavelength system. 

According to yet another preferred embodiment of the present invention, the phasor 
values of the resolution elements of the mask and the electromagnetic radiation 
25 modulator are binary, for example 0 and n. It is further preferred that the random 
distribution of the decrypting binary phasor values is uniform. 

It is also preferred that the binary values are selected so that encryption and 
decryption correspond to a binary XOR logical operation since the XOR operation 
30 provides maximum ambiguity as to which inputs generate a specific output. 

It is even more preferred that there is no spatial correlation between resolution 
elements in the encoded, non-encrypted phase pattern. 
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The reconstructed intensity pattern l(x\y') may be detected by a camera, such as a 
CCD camera, and the recorded image may be transmitted to a computer for 
recognition of the image l(x f ,y') and thus, authentication of the encrypted mask. 

5 BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 shows a 4f common path interferometer, 
Fig. 2 shows a 2f common path interferometer, 

10 

Fig. 3 shows a 1f common path interferometer. 

Fig. 4 shows (A) off-axis read-out of reflective SLM and (B) on-axis read-out of 
reflective SLM, 

15 

Fig. 5 shows schematically an example of a prescribed intensity pattern in 1D f 

Fig. 6 shows schematically the resulting phase encoding corresponding to Fig. 5, 

20 Fig. 7 illustrates a binary phase image encryption method according to the present 
invention, 

Fig. 8 illustrates a binary phase image decryption method according to the present 
invention, and 

25 

Fig. 9 illustrates alignment of an encrypted phase mask and a decrypting phase mask. 
DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

30 Fig. 1 shows a 4f common path interferometer 1 . A laser 2 emits a light beam which 
is expanded by a beam expander 3 into a plane light wave of uniform intensity and 
directs it towards an encoder 4, i.e. a combination of an encrypted mask 4a and a 
decrypting complex spatial electromagnetic radiation modulator 4b. The light beam is 
transmitted through the encoder 4 and a Fourier transforming lens 5. The encoder 4 is 
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positioned in the front focal plane of the lens 5 and a spatial filter 6 is positioned in 
the back focal plane of the lens 5 that is also the front focal plane of a lens 7. The 
Fourier transforming lenses 5, 7 need not have identical focal lengths. Different focal 
lengths lead to a magnification ratio different from one. The filter 6 phase shifts by 6 
5 and optionally attenuates (by a factor B) the zero order diffraction part 8 of the light 
modulated by the encoder 4. Optionally, the remaining diffraction part of the light 
modulated by the encoder 4 may be attenuated by a factor A. The reconstructed 
intensity pattern is generated in the back focal plane 9 of the lens 7 and a dynamic 
focusing system 10 images the reconstructed intensity pattern Kx'.y') onto a focusing 
10 plane 1 1 . The reconstructed intensity pattern l(x',y') may be detected by a camera, 
such as a CCD camera, and the recorded image may be transmitted to a computer 12 
for recognition of the image Kx'.y') and thus, authentication of the encrypted mask 
4a. 

1 5 The computer 1 2 controls the optical system. The computer 1 2 comprises interface 
means for addressing each of the resolution elements of the complex spatial 
electromagnetic radiation modulator 4b and transmitting a decrypting value to the 
addressed resolution element. Alternatively, the complex spatial electromagnetic 
radiation modulator 4b may be a fixed mask i.e. a mask manufactured with fixed 

20 modulating values of the resolution elements. 

Optionally, the phase shift 9 and attenuation factors (A, B) of the filter 6 is adjustable 
and controllable by optional phase control means of the computer 12 which may be 
further adapted to adjust the phase shift, e.g. utilising equation 18. 

25 

Fig. 2 shows a 2f common path interferometer 20. A laser 21 emits a light beam 
which is expanded by a beam expander 22 into a plane light wave of uniform intensity 
and directs it towards an encoder 23, i.e. a combination of an encrypted mask 23a 
and a complex spatial electromagnetic radiation modulator 23b and a polarisation 
30 beam splitter 24 and a quarter-wave plate 25. The polarisation beam splitter 24 and 
the quarter-wave plate 25 allows beam-splitting of light of a specific linear polarisation 
without the power loss associated with conventional beam-splitters due to splitting of 
the beam in both directions of transmission through the beam-splitter. After 
transmission through the polarisation beam splitter 24 and the quarter-wave plate 25. 
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the light beam is transmitted through a Fourier transforming lens 26 and is reflected 
from a spatial filter 27. The encoder 23 is positioned in the front focal plane of the 
lens 26 and the spatial filter 27 is positioned in the back focal plane of the lens 26. 
The filter 27 phase shifts by 9 and optionally attenuates (by a factor B) the zero order 
5 diffraction part 28 of the light modulated by the encoder 23. Optionally, the remaining 
diffraction part of the light modulated by the encoder 23 may be attenuated by a 
factor A. The reconstructed intensity pattern l(x\y'} is generated in the back focal 
plane 29 of the lens 26 and a dynamic focusing system 30 images the reconstructed 
intensity pattern l(x\y') onto a focusing plane 31. As described for the system shown 
10 in Fig. 1, the image 1(x\y') may be detected with a camera and transmitted to a 
computer 32 for processing and authentication and the system 20 may also be 
controlled by the computer 32. 

Fig. 3 shows a 1f common path interferometer 40. A laser 41 emits a light beam 
15 which is expanded by a beam expander 42 into a plane light wave of uniform intensity 
and directs it towards an encoder 43, i.e. a combination of an encrypted mask 43a 
and a complex spatial electromagnetic radiation modulator 43b. The light beam is 
transmitted through the encoder 43 and an image forming lens 44. The filter 45 phase 
shifts by 6 and optionally attenuates (by a factor B) the zero order diffraction part of 
20 the light modulated by the encoder 43. Optionally, the remaining diffraction part of 
the light modulated by the encoder 43 may be attenuated by a factor A. The 
reconstructed intensity pattern is generated in the image plane 46 of the lens 44 and a 
dynamic focusing system 47 images the reconstructed intensity pattern l(x\y') onto a 
focusing plane 48. As described for the system shown in Fig. 1, the image l(x\y') 
25 may be detected with a camera and transmitted to a computer 49 for processing and 
authentication and the system 40 may also be controlled by the computer 49. 

Fig. 4 shows details of (Fig. 4A) an off-axis read-out of an encoder 50 comprising a 
combination of a reflective mask 50a and a complex spatial electromagnetic radiation 
30 modulator 50b and of (Fig. 4B) an on-axis read-out of an encoder 51 comprising a 
combination of a reflective mask 51a and a complex spatial electromagnetic radiation 
modulator 51b with a beam splitter 52. Both configurations (Fig. 4A, Fig. 4B) may be 
utilised in the systems shown in Figs. 1-3. 
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PHASE ENCODING FOR DC PHASE FILTERING 

In the following an example of encoding an encoder and a spatial phase filter will be 
given based on a system filtering in the DC-frequency range. The exemplified system 
5 is based on a 4-f lens configuration as shown in Fig. 1 and illuminated by 

electromagnetic radiation in the visible frequency domain, hereafter simply denoted as 
fight radiation. 

Assuming that the illuminating light is monochromatic and has a substantially flat 
10 amplitude profile we obtain the following spatial amplitude distribution emitted from 
the encoder: 



15 where ct(x, y) = exp(i<t>(x, y)) represent the spatially encoded phasor values and 
AxAy is the area of the input phase modulating spatial light modulator. 

It turns out to be convenient to separate a(x, y) into two terms describing a spatially 
invariant DC-value, a , and a spatially varying AC-contribution Act(x, y) . The DC- 
20 value can be found as: 




(1) 




(2) 



Subsequently the AC-term is expressed by: 



25 



Act(x, y) = expfitKx, y)) - — ^T" JJ exp(i(Kx, y))dxdy (3) 

AxAy ^ y 
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The separation of ct(x, y) into a spatially invariant DC-term and a spatially varying 
AC-term is an important point and will be used throughout the remaining part of this 
example, especially in the description of the spatial filtering procedure. 

5 The spatial filter utilised in this example is chosen as a circular phase contrast filter 
{different transverse shapes can also be used) centred around origo in the spatial 

frequency domain, denoted by co-ordinates (f x , f y ): 



T (f r ) = 1 + (exp(ie) - l)c 




10 

where f r = ^Jf* + f y 2 denotes radial spatial frequency and Af r describes the size 
of the circular (circ) phase filter. 



In the spatial frequency domain (the filtering plane) the Fourier transformation ( 3) of 
1 5 the spatially modulated light radiation from the encoder is present. The filtering 
operation on the Fourier transformed light radiation performed by the spatial phase 
contrast filter can be expressed as a simple point-by-point multiplication procedure. 
Subsequently the spatially filtered light is inverse Fourier transformed ( X 1 ) by the 
second Fourier lens (Fourier transformation and reflected output co-ordinates) and the 
20 resulting spatial amplitude distribution in the image plane (with co-ordinates (x' f y') ) 
can accordingly be written as: 



f 

o(x', y') = a(x>, y') + (exp(i9) - l)^ 1 



3(a(x, y))circ^j 
= a(x', y') + a(exp(i0) - l)rect(^ , -^j (5) 
= [exp(i<Kx\ y')) + a(exp(ie) - l)]rect^ , 
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Within the illumination-region, (x f , y') <= W , outlined by rect ^— , ^— , 

V Ax Ay ) 



one 



obtains: 



|o(x', yf = 1 + sin (f sin (|) " sin (*i " ^ X '' y,) + I)] (6) 

5 

Requiring that |oC^, y'^j[ — 0 corresponding to complete darkness as the lowest 
intensity level in regions (x' ot y' a ) e H' 0 implies: 



1 + 

10 



fl sin(|||a| sin[|) - sinj^ - 4> 0 + |j] = 0 (7) 



where the abbreviation <{> 0 = §(x r 0 , Yo) has been used - 
The solutions to Eq. (7) are given by: 



sinf*- -*. + |)± Jsin 2 (^ - * e + |) - 
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1 



2 S1I 



The requirement o < |a| < 1 implies that: 



(8) 



G = n - 2((|>- - <|> 0 ) + p2rc, p = 0,114:2, . . . 

20 

leading to 



20 



where the +sign is for 0 -values in the interval: 




(11) 



5 

and the -sign is for 9 -values: 

e e + p«»27i (12) 

10 The corresponding interval for - <J> 0 ) is: 

(♦.--♦.) e ]i'-f[ (13) 

Inserting the expression for |ot| , one obtains the simple intensity expression: 

15 

|o(x', yf = 2^1 T s±t^- - <Kx', Y') + I)] d4) 

where 

JJ exp(i<|>(x, y))dxdy = AxAy|oc| exp(i<|>-) (15) 

AxAy 

20 The phase-only transformations imply that energy is conserved: 

JJ |o(x', y'fdx'dy' = JJ |a(x, y)fdxdy = AxAy (16) 

25 A special case: 
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The most convenient choice for a is: a = — (implying that 6 = n + p even 27i ), so 

2 

that the output intensity can be described as: 

|o(x\ y')f = 2[l - cos(<Kx', y'))] "7) 

5 

In this case the phases-intensity mapping is described by the intervals 
[0; it] -> [0;4] . 

By setting a = — one obtains the following requirement to the phase function 
y) : 

10 

AxAy 



JJ cos($(x, y))dxdy = 



AxAy " (18) 

JJ sin(<|>(x, y))dxdy = 0 

AxAy 



Inserting the expression for |o(x', y')f in Eq. (16) yields: 
15 2 JJ [l - cosfo(x\ y'))fc*'dy' = AxAy (19) 

AxAy 

in accordance with the first of the integral expressions in Eq. (18). 



20 

Encoding procedure: 

• A given intensity distribution (image) |o(x', y')| 2 is desired at the output side of 
the optical set-up. 

25 



15 
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Pixellation of the image, that is generally represented in the greyscale range: 
[0; gmax] , provides the relation: 

f f |o(x', y'fdx'dy' = AxAy => £ H 1 ' j)f = P lj W) ' 



The histogram for the desired image |o(i, j)f is adjusted (adj.) within the 
greyscale range [0; gmax], so that the previous point is fulfilled: 

fast -> MMU- 



The phase values can now be calculated as: 



10 <|)(/,y) = arccos 



gmax 



J 



As before pixellation provides the relation: 
£ sin((Ki, j)) = °- 



• The previous point can now be fulfilled by complex conjugating half the input 
pixels having the same phase value in the phase histogram. 

• The phase conjugate phase flipping provides a valuable tool (an extra degree of 
20 freedom) for manipulating the spatial frequency content in order to optimise the 

separation of low and high frequency terms at the filter plane. 

• The scheme is robust to constant phase errors across the input spatial phase 
modulator, since Eq. (14) is a function of the difference: <^ o - 4>(i, j) * on| V. 

25 Furthermore, small variations in the individual pixel phase values do not introduce 

any detrimental effects because the average value a , is a result of a very large 
phasor sum. 
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• if the desired intensity distribution is too small to include all energy, that is, the 
histogram is scaled to maximum and the left hand side of Eq. (16) is still smaller 
than the right hand side, then the input phase object can be scaled until Eq. (16) 
is fulfilled. In order to obtain a scale invariant output intensity level a dynamic 
5 focusing system is needed. Similarly, intensity invariance can be obtained by 

controlling the radiated power from the light source. Alternatively, one can ignore 
the residual background illumination and obtain intensity levels with a gain factor 

of 9" (background constant equal to 1") for narrow generally shaped line 
structures (e.g. Eq. (6)). 

10 

Example 1: 

A very simple example illustrating the individual steps in the above procedure will be 
given below. To simplify the example it will be considered in one dimension only. The 
starting point for encoding the encoder in this example is based on the following 
15 parameters: 



a = 0.5 

6 - * (20) 

gmax — 4 



Consider the pixellated 3-step function shown in Fig. 5 to be reconstructed in the 
20 image plane as an intensity distribution. From the above choices of parameters one 
obtains the simple relation between phase values in the encoder and the image 
intensity values: 



\o(if = 2[l - cos(*(i))] (21) 

25 

To proceed from here it necessary to calculate the accumulated intensity ^ |°(^)] 2 ' n 

i 

the image to be reconstructed. The accumulated intensity is easily calculated from an 
image histogram where the x-axis represents greylevel value and the y-axis represents 
the amount of pixels in the image at a given greylevel value. By use of a histogram 
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£ \o(if is simply found as the weighted sum of all greylevel values (x-axis) 
i 

multiplied by their pixel counting (y-axis). This describes, so to speak, the "weight" of 
the image. In this simple example histogram calculations are not needed since we only 
have 3 greylevels with well-defined separations. 

The value for the accumulated intensity has to obey the equality: 



£ |o(i)| 2 = ^# pi* (Ajt) =# pi*^, = 1 



= 14 



(22) 



10 From Fig. 5 we obtain: 

£ |o(i)| 2 = 4pixels • 0 + 4pixels • (OSmax) + 6pixels • max = Smax 

i 

(23) 

15 So that the value for max can be estimated to be: 



max — — 
4 



(24) 



The corresponding adjusted intensity levels, |o(ij|* dj , are therefore: 7/4, 7/8 and 0. 
20 These values can now be utilised to calculate the phase values of the encoder from 
the relation: 



<|>(i) = arccos 



2 \ 



1 ~ 



gmax 



= arccos 



1 - 



H4 



2 \ 

adj 



(25) 



25 where from we obtain the three phase values: 1.45 rad. 0.97 rad. and 0 rad. 

The last step needed in order to encode the encoder is that the following equality is 
fulfilled: 



25 

£ sin/>(i)) = 0 (26) 
i 

Since we have the choice to use complex conjugate phasor values (two phasors giving 
the same intensity level) many approaches can be taken from here. A simple approach 
5 is to flip every second phasor with its complex conjugate value as shown in Fig. 6. 
The final phase values used in the encoder are accordingly: ± 1 .45 rad. ± 0.97 rad. 
and 0 rad. 

As the last step we can check whether the criteria: a = 1 / 2, is actually fulfilled 
1 0 with the chosen phasor encoding: 



a = — (4 exp(iO) + 2 exp(i0.97) + 3 exp(il.45) + 2 exp(-i0.97) + 3 exp(-il-45)) s 1 / 2 
14 v 

(27) 



1 5 GENERAL PHASE CORRECTION PROCEDURE INTEGRATED WITH THE PHASE 
ENCODING 

In Eq. (6) we obtained an analytic relation between the phase values in the encoder 
and the resulting intensity distribution, within the region (x", y') e 91' : 

20 

|o(x', yf * 1 + «R *Kf 1^ ""(I) " Sin (*«- " ^ + 2)] (28> 

The analysis leading to the above relation was based on the assumption that |~c| is a 
constant value within the 9T -domain. In other words, the following approximation 
25 was applied: 



However, for certain spatial filter parameters the left-hand side of this expression wi 
not be a space invariant constant value throughout the whole 9T -domain but will 



26 



instead manifest slowly variations/oscillations. This will introduce small errors in the 
final superposition between the phase filtered DC-value and the direct propagated AC- 
signal. In order to circumvent this problem a technique is needed that can counteract 
the distortions by use of phase-only encoding in the components already present in 
5 the system. In what follows a procedure for integrating predistortion that counteracts 
the above mentioned distortions will be described that is purely based on modifying 
the phasor values in the encoder at the input side of the system. The method can also 
counteract other types of distortions inherent in a practical implementation of the 
system. Furthermore, the method can be applied in systems filtering at other spatial 
10 frequencies than DC. 

Procedure: 

When encoding the input phase function it is helpful to have a 'reverse' equation, 
expressing the input phase distribution as a function of an adjusted (electronic) image 
1 5 grey-level distribution, I Blm , addressing the input spatial light modulator: 

S 1 + ^(x',y)|sin[|Jp(x' /y ')|sin(|) - sin^x-, y') - ♦(x\y) + |)] 



gmax 

(30) 

where it has been taken into account that a(x', y') is not considered as a constant 
20 but manifests a smooth oscillating behaviour within the optical image domain. The 
maximum value of J Blm is denoted gmax . 

Now, one can derive a formula for the 'grey-level correction' Ar slm (x', y') that one 
needs to apply in order to encode a phase function that compensates for the spatial 
25 variation of the average phase value a(x\ y') : 

«*..■<*'•>'•> = i + ^(x' ( y')| S in(||F(x', y')| sing) - sin(* s (*\ y') - ♦<*•. yO + f )] 



gmax ^ 
/ A r 2(l.,.(x',y') + AJ.^x'.yOn 
t (x« f y') = acos^i — J 



(31) 
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where the second relation has been derived from the first by setting a =-j and 
0 = 7i . 

By inserting the second relation in the first expression one gets: 



5 

AJ sl „(x\ Y') = 



(32) 



This formula is however not directly useful because it is related to the histogram 
1 0 adjusted grey-level distribution denoted by I slm - 

One needs a formula that relates the above correction term to the 'original' input grey- 
level distribution J(x, y) that has not been modified by histogram adjustments. This 
is important since the effect of the grey-level corrections also have to be incorporated 
15 in the procedure of histogram adjustments. 

The histogram scaling gives: 

J(X/ y) = i slm {x, y) (33) 

20 



where J max and I almtmBX are the maximum grey-level values occurring in the original 
and the adjusted electronic grey-level distributions respectively. 

Similarly, one can apply this relation to the intensity correction term A I glm and obtain: 

25 

T(x, y) = J(x, y) + A J(x, y) = Jmax (x alJ .(x f y) + M elm (x, y)) 

sim.max 

(34) 

resulting in: 
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f(x - * ■ ^ [ iix - y) - w ^ - »T) l36 ' 

In order to have 'enough dynamic range' in grey-levels for the correction term one can 
5 derive an inequality from the above relation by using the fact that J max < gmax: 



loin ^ ^slnMnax ^ ' / 



or 



2ja ra in| • gmax 



(37) 



10 

Since the first term is the dominating term in the expression for the intensity 
correction it will in practice be sufficient just to have the much simpler corrections: 



15 



y ^ _ y) 
2fx(x, y)| 
* 2 |rl mln | • gmax 

BINARY PHASE IMAGE ENCRYPTION AND DECRYPTION METHOD 



(38) 



Fig. 7 illustrates a binary phase image encryption method and Fig. 8 illustrates a 
binary phase image decryption method according to the present invention. In the 
20 illustrated example the modulating phase values are 0 or n. In Figs. 7 and 8 bright 
resolution elements 100 have the phase shift value 0 and the dark resolution elements 
101 have the phase shift value jr. 



The encryption method comprises the steps of 

25 
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1) calculating, according to the mathematical method described above and in WO 
96/34307, the modulating phase values (0,7i) of the resolution elements forming a 
binary encoded phase pattern 102 that is imaged onto the intensity pattern 103, 
e.g. by one of the systems illustrated in Figs. 1-3, 

5 

2) generating a uniform random distribution 104 of binary phase values (0,7r) r 



3) forming an encrypted phase mask 105, 4a, 23a, 43a, 50a, 51a by adding, for 
each resolution element 106 of the encrypted mask 105, the calculated phase 

10 value 107 (0 or 71) and the corresponding generated random phase value 108 {0 or 
n) of the resolution element, whereby a uniform random distribution of binary 
phase values (0,7i) is generated in the encrypted phase mask 105 encrypting the 
image 103, and 

1 5 the decryption method comprises the step of reconstructing the image 103 by 

4) adding the binary phase. shift value (0,ti) of each resolution element of the 
encrypted phase mask 105 with a binary phase shift value (0,ti) of a corresponding 
resolution element of a decrypting phase mask 109, 4b, 23b, 43b, 50b, 51b with 

20 resolution elements having the respective calculated phase values 104, whereby 
the original encoded phase pattern 102 is reconstructed. 

The decrypting step 4) may be performed by aligning the encrypted phase mask 105 
and the decrypting phase mask 109 in one of the optical systems illustrated in Figs. 1- 
25 3, whereby the original image 103 is reconstructed, e.g. for detection by a camera 
and recognition by a computer. 

It should be noted that in the above example, encrypting and decrypting phase values 
of corresponding resolution elements are identical. 

30 

Utilisation of 0 and n as the binary phase values provides a particular robust 
encryption and decryption approach with a low sensitivity to perturbations of the 





"Modtaget PD 
15 MRS. 
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masks because of the relatively large phase difference between the binary values of 



In the common path interferometer illustrated in Fig. 9, the decrypting mask 1 20b is a 
5 spatial light modulator with a modulating region that is larger than the encrypted 
phase mask 120a. Thus, the encrypted phase mask 120a and the decrypting mask 
120b can be aligned electronically by appropriate control of the spatial light modulator 
120b. For example, the encrypted mask may contain regions with non-encrypted 
patterns that are imaged onto the CCD detector array 121 for the purpose of 
10 alignment. The positions of the imaged alignment patterns are detected and the spatial 
light modulator 120b is subsequently controlled accordingly whereby the decrypting 
resolution elements of the light modulator 120b are aligned with the resolution 
elements of the encrypted phase mask 120a. 

15 

CLAIMS 

1 . A method of decryption of an encrypted image having a non-encrypted image 
intensity pattern l(x\y') and 

20 

Encoded into a mask having a plurality of mask resolution elements (x m ,y m ) with an 
encoded phase value 4>(x m ,yJ and an encoded amplitude value a(x mf y m ), and 

Encrypted by addition of an encrypted phase value 4> c (x m ,Y m ) to the encoded phase 
25 values <|>{x m ,yj and by multiplication of an encrypted amplitude value a c (x m ,yj with 
the encoded amplitude value a(x m ,yj, 

Each mask resolution element (x m ,y m ) modulating the phase and the amplitude of 
electromagnetic radiation incident upon it with the complex value 



the mask. 



30 a(x, 



Vm)a c (x ( 




The method comprising the steps of 
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Radiating electromagnetic radiation towards the mask, 

Inserting into the path of the electromagnetic radiation a complex spatial 
electromagnetic radiation modulator comprising modulator resolution elements (x d ,y d >, 
5 each modulator resolution element (x d ,y d ) modulating the phase and the amplitude of 
electromagnetic radiation incident upon it with a predetermined complex value 
a d (x d ,y d )e i ♦ d,xd ■ vd, , the decrypting phase value <|> d (x d ,y d ) and the decrypting amplitude 
value a d (x d ,y d ), respectively, of a modulator resolution element <x d ,y d ) being 
substantially equal to -<Mx m ,yJ and a c 1 (x m ,yj # respectively, of a corresponding mask 
10 resolution element (x m ,y m ), and 

Imaging the mask and the electromagnetic radiation modulator onto the image having 
the image intensity pattern l(x\y"). ' 

15 2. A method according to claim 1, wherein the step of imaging comprises imaging 
with a common path interferometer. 

3. A method according to claim 1 or 2, wherein the step of imaging comprises phase 
contrast imaging. 

20 

4. A method according to claim 2 or 3, further comprising the steps of 

Fourier or Fresnel transforming electromagnetic radiation modulated by the mask and 
the complex spatial electromagnetic radiation modulator, 

25 

Filtering the Fourier or Fresnel transformed electromagnetic radiation by 

In a region of spatial frequencies comprising DC in the Fourier or Fresnel plane, 



30 



Phase shifting with a predetermined phase shift value 6 the modulated 
electromagnetic radiation in relation to the remaining part of the 
electromagnetic radiation, and 
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Multiplying the amplitude of the modulated electromagnetic radiation 



with a constant B, and 



5 



In a region of remaining spatial frequencies in the Fourier or Fresnel plane, 
multiplying the amplitude of the modulated electromagnetic radiation with a 



constant A, 

Forming the intensity pattern by Fourier or Fresnel transforming, respectively, the 
phase shifted Fourier or Fresnel transformed modulated electromagnetic radiation, 
10 whereby each resolution element (x m ,y m ) of the mask is imaged on a corresponding 
resolution element Ix*^') of the image, 

The filtering parameters A, B, 0 substantially fulfilling the equation 



a being the average of the complex phasors a{x m/ y m )e i ** Jtrn ym \ 

5. A method according to claim 4, wherein the filtering parameters A and B 

20 substantially fulfil that A = 1 and B = 1 , and wherein the absorption of the mask is 
substantially uniform. 

6. A method according to claim 5, wherein the phase shift value 9 substantially fulfils 
the equation 



15 



l(x\y') = A^atx'^^e^^ + BA' 1 a (e iB -1)| 2 



25 




1 



7. A method according to claim 6, wherein the phase shift 9 is substantially equal to 
30 7i. 
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8. A method according to any of claims 3-7, further comprising the steps of 

Moving the DC-part of the electromagnetic radiation to a second part of the Fourier or 
Fresnel plane, and 

5 

Phase shifting the Fourier or Fresnel transformed modulated electromagnetic radiation 
at the second part of the Fourier or Fresnel plane by G in relation to the remaining part 
of the electromagnetic radiation. 

10 9. A method according to claim 8, wherein the step of moving the DC-part of the 
electromagnetic radiation comprises utilisation of an optical component, such as a 
grating, a prism, etc, with an appropriate carrier frequency. 

10. A method according to any of claims 3-9, further comprising the step of phase 

1 5 shifting at selected spatial frequencies constituting a region that is shaped to match 
the spatial frequency con-tent of the phasors e* ,xm,Yml . 

1 1 . A method according to any of claims 3-10, wherein the step of filtering comprises 
utilisation of a spatial light modulator. 

20 

1 2. A method according to any of claims 3-1 1 , further comprising the step of 
encoding the optical function of an output lens into the filter. 

13. A method according to any of claims 3-12, wherein the step of radiating 

25 electromagnetic radiation comprises radiation of electromagnetic radiation of different 
wavelengths corresponding to three different colours, such as red, green and blue, for 
generation of intensity patterns of arbitrary colours. 

14. A method of encryption of an image having an intensity pattern l(x\y') to be 
30 decrypted according to any of claims 1-15, comprising the steps of 



Pixellating the intensity pattern l(x\y') in accordance with the disposition of resolution 
elements (x m ,y m ) of a mask, 
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Encoding the mask with an encoded phase value (t>(x m ,yj and an encoded amplitude 
value a{x m ,y m ), and 

Encrypting by addition of an encrypted phase value $ c (x m ,y m ) to the encoded phase 
5 values 4>(x m ,y m > and by multiplication of an encrypted amplitude value a c (x m ,yj with 
the encoded amplitude value a(x m ,y m ), 

Each mask resolution element (x m ,y m ) modulating the phase and the amplitude of 
electromagnetic radiation incident upon it with the complex value 
10 a(x m( y m )a c (x mf y m )e^^ + ^'. 

1 5. A method according to claim 14, further comprising the step of calculating the 
complex phasor values a{x m ,yje* ( xm ' ym) of the mask in accordance with 

15 Hx\y') = A^atx^'Je*^ + BA" 1 a {e^-DI 2 

For selected phase shift values 8, a being the average of the complex phasors 

a/v w \ 0 '*l xm.ym) 
aix m ,y m /e # 

20 Selecting, for each resolution element, one of two phasor values which represent a 
particular grey level. 

16. A method according to claim 15, wherein the filtering parameters A and B 
substantially fulfil that A = 1 and B = 1 , and wherein the absorption of the mask is 

25 substantially uniform. 

17, A method according to claim 16, wherein the step of calculating the phasor 
values comprises 

Setting the reconstructed intensity of at least one resolution element (x 0 \ y 0 *) of the 
30 intensity pattern to zero, and 



Calculating the phasor values e**-* of the mask in accordance with 
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N = ] -k- 

2|sin|| 



/(*' ,y ) = 2^1 + sin(<|^ -Mx'>y' ) +|)] 



For selected phase shift values <|>, being the phase of a . 

10 18. A method according to claim 17, further comprising the step of selecting the 
phase shift 9 = n, selecting |a| = Vi . and calculating the phasor values e* (xy > of the 
encoder in accordance with 

Kx'.y') = 2(1 -cost Wx-.y'))) 

15 

jjsin(b(x m ,yjdx m dy ni = 0 . 

mask 



19. A method according to any of claims 14-18, further comprising the step of 

20 encoding the function of an optical component, such as a grating, a prism, etc, with 
an appropriate carrier frequency, into the mask. 

20. A method according to any of claims 14-19, further comprising the step of 
adjusting the modulus of the Fourier transform of the complex phasors a(x m ,y m )e* ,xm ' VfT,, 

25 at specific spatial frequencies in order to control the range of intensity levels of the 
reconstructed intensity pattern. 
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21 . A method according to claim 20, wherein the step of adjusting the modulus of the 
Fourier transform of the complex phasors a(x m ,yJe* ,XTnyffl> at specific spatial 
frequencies comprises at least one of the following measures: 

5 a) Adjusting the individual complex phasors a<x m ,yje* ,xm ' Ym) of the resolution 

elements of the mask maintaining prescribed relative intensity levels between 
intensities of resolution elements of the intensity pattern, 

b) Adjusting the individual complex phasors a<x m ,yje* ( xm * vm> of the resolution 
10 elements of the mask by histogram techniques, 

c) Spatially scaling the complex phasor a(x m ,y m )e i * txm vni) pattern of the mask, and 

d) Utilising half tone coding techniques. 

15 

22. A method according to any of claims 14-21, wherein each complex phasor 
a(x m ,yje* txm ' vm> Q f t ^ e mask is selected from a set of two determined phasors with 
complementary complex phasor values a{x m , Ym )e i * Mxm vm) and atx^y^* 2 ^'^ in such a 
way that a specific spatial frequency distribution of the intensity of the 

20 electromagnetic radiation in the Fourier or Fresnel plane is attained. 

23. A method according to claim 22, wherein the phase <|><x,y) of complex phasors 
a(x m ,yje i * txm * vm) of adjacent resolution elements alternates between the two possible 
complementary complex phasor values ^(x m ,yJe* Mxm • vm, and a{x m ,yJe* 2(xm < vm ^ 

25 

24. A method according to claim 22 or 23, wherein the complex phasors 
a(x m ,y m )e i * 1{xm * vm) and a(x m ,yJe , * 2Urn ' vm » are complex conjugated. 

25. A method according to any of claims 14-24, further comprising the step of 

30 encoding the optical function of a Fourier-transforming lens into the complex phasors 
a(x m ,yje* (xm ' vm) of the mask. 
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26. A method according to any of claims 14-25, wherein each of the encrypted phase 
values $ c (x m ,y m ) is substantially equal to a value selected from a set consisting of two 
phase values. 



5 27. A method according to claims 26, wherein each of the encrypted phase values 
(() c {x m ,y m ) is substantially equal to a value selected from the set consisting of 0 and rc. 

28. A decryption system for decrypting an encrypted image having a non-encrypted 
image intensity pattern l(x\y ? ) that has been 

10 

Encoded into a mask having a plurality of mask resolution elements (x m# y m ) with an 
encoded phase value (J>(x m ,y m ) and an encoded amplitude value a(x m ,y m ), and 

Encrypted by addition of an encrypted phase value <|> c (x m ,y m ) to the encoded phase 
15 values <|>(x m ,y m ) and by multiplication of an encrypted amplitude value a c (x m ,y m ) with 
the encoded amplitude value a(x m ,y m ), 

Each mask resolution element (x m ,y m ) modulating the phase and the amplitude of 
electromagnetic radiation incident upon it with the complex value 
20 a(x m/ y m )a c (x mf y m )e rtum ' vm)+i * cumvm, f 

The system comprising 

A source of electromagnetic radiation for emission of electromagnetic radiation, 

25 

Positioned in the path of the electromagnetic radiation a complex spatial 
electromagnetic radiation modulator comprising modulator resolution elements (x d ,y d ), 
each modulator resolution element (x d ,y d ) modulating the phase and the amplitude of 
electromagnetic radiation incident upon it with a predetermined complex value 
30 a d (x d ,y d )e* dud vd) , the decrypting phase value <f> d (x d ,y d ) and the decrypting amplitude 
value a d (x d ,y d ), respectively, of a modulator resolution element (x d ,y d ) being 
substantially equal to -<{> c (x m ,yJ and a c '*{x m ,y m ), respectively, of a corresponding mask 
resolution element (x m ,y m ), and 





38 

An imaging system for imaging the mask and the electromagnetic radiation modulator 
onto the image having the image intensity pattern l{x\y f ). 

29. A system according to claim 28, wherein the imaging system comprises a 
5 common path interferometer. 

30. A system according to claim 28 or 29, wherein the imaging system comprises a 
phase contrast imaging system. 

10 31 . A system according to claim 29 or 30, further comprising 

Means for Fourier or Fresnel transforming the electromagnetic radiation modulated by 
the mask and the complex spatial electromagnetic radiation modulator and being 
positioned on a propagation axis of the modulated radiation. 



15 

A spatial filter for filtering the Fourier or Fresnel transformed electromagnetic radiation 



by 



In a region of spatial frequencies comprising DC in the Fourier or Fresnel plane, 



20 



Phase shifting with a predetermined phase shift value 9 the modulated 
electromagnetic radiation in relation to the remaining part of the 
electromagnetic radiation, and 



25 



Multiplying the amplitude of the modulated electromagnetic radiation 
with a constant B, and 



In a region of remaining spatial frequencies in the Fourier or Fresnel plane, 
multiplying the amplitude of the modulated electromagnetic radiation with a 



30 



constant A, 



Means for forming the intensity pattern by Fourier or Fresnel transforming, 
respectively, the phase shifted Fourier or Fresnel transformed modulated 




39 

electromagnetic radiation, whereby each resolution element (x m ,y m ) of the mask is 
imaged on a corresponding resolution element (x\y T ) of the image, 

The filtering parameters A, B, 0 substantially fulfilling the equation 

5 

Kx',y') = A^aU'^y'Je^' 1 + BA 1 a (e i9 -1)| 2 

For selected phase shift values 9, a being the average of the complex phasors 
a(x m ,yje* (xmvm \ 

10 

32. A system according to claim 31 , wherein the filtering parameters A and B 
substantially fulfil that A = 1 and B = 1 , and for each (x m ,y m ) of the mask: a{x m ,yj = 1 . 



33. A system according to claim 32, wherein the phase shift value 8 substantially 
1 5 fulfils the equation 



|a| = - 



2|sin ■ 



34. A system according to claim 33, wherein the phase shift 0 is substantially equal 
20 to 7i. 

35. A system according to any of claims 30-33, further comprising 

Means for moving the region of spatial frequencies comprising DC to a second part of 
25 the Fourier or Fresnel plane, and wherein 

The spatial filter is adapted to phase shift the transformed modulated electromagnetic 
radiation at the second part of the Fourier or Fresnel plane by 0 in relation to the 
remaining part of the electromagnetic radiation. 



30 




40 

36. A system according to claim 35, wherein the means for moving the region of 
spatial frequencies comprising DC to a second part of the Fourier or Fresnel plane 
comprises an optical component, such as a grating, a prism, etc, with an appropriate 
carrier frequency. 

5 

37. A system according to any of claims 31-36, wherein the spatial filter comprises a 
spatial light modulator. 

38. A system according to any of claims 31-37, wherein the spatial filter is adapted to 
10 perform the optical function of an output lens by appropriate encoding of the spatial 

filter. 

39. A system according to any of claims 28-38, wherein the source of 
electromagnetic radiation is adapted to radiate electromagnetic radiation of different 

1 5 wavelengths corresponding to three different colours, such as red, green and blue, for 
generation of intensity patterns of arbitrary colours. 

40. A system according to any of claims 31-39, further comprising a first and a 
second Fourier transforming lens, the mask being positioned in the front focal plane of 

20 the first lens, the spatial filter being positioned at the back focal plane of the first lens, 
and the second lens being positioned so that its front focal plane is positioned at the 
position of the back focal plane of the first lens. 

41 . A system according to any of claims 31-40, further comprising one Fourier 

25 transforming lens, the spatial filter being positioned at the back focal plane of the lens. 

42. A system according to any of claims 31-41, further comprising one imaging lens, 
the spatial filter being positioned in the back focal plane of the lens. 

30 43. A system according to any of claims 31-42, further comprising a polarising beam 
splitter and a quarter wave plate and/or a phase filter reflecting electromagnetic 
radiation incident upon it. 



41 

44. A system according to any of claims 31-43, wherein the spatial filter changes the 
phase of the radiation in the region of spatial frequencies comprising DC and leaves 
the phase of the remaining part of the radiation unchanged. 

5 45. A system according to any of claims 31-43, wherein the spatial filter do not 
change the phase of the radiation in the region of spatial frequencies comprising DC 
and changes the phase of the remaining part of the radiation. 

46. A system according to any of claims 31-43, wherein the spatial filter blocks the 
10 radiation at the region of spatial frequencies comprising DC and leaves the remaining 

part of the radiation unchanged. 

47. A system according to any of claims 31-46, wherein the source of 
electromagnetic radiation is a Laser. 

15 

48. A system according to any of claims 28-47, wherein each of the encrypted phase 
values <Mx m ,y m ) ' s substantially equal to a value selected from a set consisting of two 
phase values. 

20 49. A method according to claims 48, wherein each of the encrypted phase values 
<|> c (x m ,y m ) is substantially equal to a value selected from the set consisting of 0 and n. 
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a) encoded phase pattern b) encrypting phase pattern c) encrypted phase mask 

Encryption Process for a binary phase image. The phase value of the pixels is either 0 (shown as light 
grey areas) or n (dark grey). The resulting phase mask (c) is a uniform random distribution of binary phase 
values and is generated by a multiplication of patterns (a) and (b) within the computer. 



Fig. 7 
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a) encrypted phase mask b) decrypting phase mask c) decrypted phase pattern 



d) intensity readout 



The decryption process. This process is the reverse of that shown in Fig. 1, but is carried out 
optically, the high contrast image (d) is obtained using a n phase contrast filter in the decryption optics. 



Fig. 8 
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